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Abstract. We consider a pointwise stabilization problem for a coupled wave 
and plate equations. We prove under rather general assumptions, that such 
systems can stabilized so as to have arbitrarily high decay rates and are exactly 
controllable. We propose a numerical approximation of the model and we study 
numerically the construction of the feedbak law leading to exponential decay 
with arbtrarily large rate. 



1. Introduction 



Let £,77 £ (0, 7r) are given points, fix four real numbers A, £?, C, D and con- 
sider the coupled string-beam system, more precisely we have the following partial 
differential equations with pointwise dissipation: 



(1.1) 



Vi,tt ~ Vi.xx + Ayi + Cy 2 = vi(t)S ( 
V2,tt + y2.xxxx + By 1 + Dy 2 = v 2 (t)S ri 
yi (t,0)= yi (t,7r)=0 
y 2 (*,0)=j/ 2 (t,7r)=0 
V2,xx{t,0) = y2,xx(t,Tr) = 
?/i(0, x) = yw(x) and j/i, t (0, x 
.2/2(0,2;) = y 2 o(x) and y 2 ,t(0,x 



in E x (0,tt), 
in E x (0,tt), 
for t E E, 
for t E E, 
for t E E, 
yii (x) for x € (0,7r), 
2/21 {x) for x € (Q,7r) 



(K), S, 



and 6 V denotes the Dirac 



where v\(t),v 2 (t) are the control functions in L 
mass at some given points £ and rj G (0, it). 

The coupled structual model has been of great interest in recent year; for details 
about the physical motivation for the model see [TH] , [E] and the references therein. 
Mathematical analysis of coupled partial differential equations is detailed in [T5] . 
|20j . The question of controllability and stabilization for such models has been 
widely treated in a series of relevant works [J , [2] , g] , 0] , Q3], [T5], [IS], [2J]- Many 
works were devoted to the construction of explicit feedback laws and to the proof 
of exponential decay by different methods; see, e.g., [5D], [13], [IS]. It is known that 
this type of feedback does not yield arbitrarily large decay rates. It was pointed out 
earlier by Haraux and Jaffard , [10] , [II] that the observability and controllability 
properties depend heavily on the location of the observation of control point. For 
the stabilization another difficulty appears because the suitable function spaces, as 
we will show, are not Sobolev spaces. 

In this paper we apply another approach for the stabilization of the coupled 
string-beam system ( |1.1[ ), which similar in sprit to the HUM. This method, devel- 
oped by Komornik [15] is as general as the former one; however, it provides stronger 
results with simpler proofs and it's the first time that we apply this method to prove 
the exponential stability of the coupled sting- beam systems with pointwise con- 
trol. The main result of this paper is to introduce functions spaces depending on 
the arithmetical properties of the stabilization point and to give an estimate on 
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exponential decay that is valid for regular initial data, as a result we will construct 
pointwise feedbacks leading to arbitrarily large prescribed decay rates. Numerical 
tests and explicit construction of the feedback are presented. 

The method used is based on a regularity results combined with an observability 
inequality for the corresponding undamped problem. See |15j . 

The paper is organised as follows. The statement and the proof of the main 
results are given in the sections 2 and 3 respectively. The last section is devoted 
to the numerical approximation of the coupled string-beam system and the explicit 
construction of the feedback law. 

2. Statement of the main result 

In order to formulate our result, we assume that !;/it and n/n are irrational, so 
that sin fc£ and sin krj, don't vanish for any k = 1,2,..., we denote by Z the linear 
hull of the functions Wk{x) :— -\/2/tt sin kx, k — 1, 2, ... , and we denote by and 
(Da)' for every a <E K and /? <E {£, n} the Hilbert spaces obtained by completing 
Z with respect to norms given by the following formulae: 

\\^2a k w k \:=Y J k 2a siu 1 (kp)\a k \ 2 , 

\zZ akWk f (Da y : = H fc_2Q sin- 2 (fc/3)|a fc | 2 . 
If we identify L 2 (0, ir) with its dual and take into account that 



\\J2 akWk ^ 2 (o^) ■= 



ak\ 2 , 



then (D^Y is the dual space of D^ . 

Fix £,77 € (0, 7r) such that £/ir and v/t: are irrational and introduce the Hilbert 
space 

H e , n := (D° x D- 1 )' x (£>° x D~ 2 )'. 



The problem ( 1.1 ) is well posed in the Hilbert space in the following sense: 
Proposition 2.1. for any given initial and final data 

(2/10,2/11,2/20,2/21) e H^ ri 

and 

v 1 ,v 2 EL 2 (0,T;H^ 1 ) 
The system has a unique weak solution satisfying 

(2/1,2/1,4,2/2,2/2,4) S C([0,T];^,^), 

and the linear mapping 

(2/10, 2/u, 2/20, 2/21, Vi,v 2 ) (2/i,2/i,t,2/2,2/2,i) 
is continuous with these topologies. 

We shall study the controllability of the system. 

Definition. Fix £, n € (0,7r) such that £,/tt and n/ir are irrational, the system 



(1.1) is exactly controllable if for any given initial and final data 

(2/10,2/11,2/20,2/21) e "H^ r) 

and 

(^10, ^11, z 20i z 2l) & H$,r] 

there exist control functions 

V!,V2 € L 2 (0,T;U^) 
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such that the corresponding solution of satisfies the final condition 

(yi,yu,y2,y2t)(T) = (z 1Q ,z 11 ,z 2 o,z 2 i). 



Theorem 2.2. If T > 2tt, then the system (1.1) is exactly controllable for almost 
all choices of (A, B, C, D) G R. 

Finally, we are looking for stabilizing feedbck laws of the form 

(2.1) {vi,v 2 ){t) := {(Pmt + Qiyi)(t,0>( p 2V2t + Q2y2){t,v)), 

leading to arbitrarily high decay rates. 

Theorem 2.3. Fix £ and rj G (0, ir) such that £/ir and rj/iv are irrational. For 
almost all choices of {A, B, C, D) G R 4 and for every positive number u> there exist 
two linear operators 



(Pi,Qi,P2,Q 2 ) : n Ln -+ DJ 1 x D- 2 , 



and a positive constant M such that the problem ( 1.1 ) is well posed in and its 
solutions satisfy the inequality 

(2-2) ||(j/i, j/it, J/2, J/2t)||H e ,„ < Me~ ut \\(y 10 , y n , y 2a , J/2i)||« Si „ 

for all (j/io,J/n, J/20, J/21 ) € Ht, v and t > 



Remark. It follows from some results of Komornik and Loreti that the system ( 1.1 ) 
can not be exactly controllable for some exceptional choices of the parameters A, 
B, C, D: see |16j and |17| for explicit counter examples concerning an equivalent 
observability problem. 

The proofs are based on the study of the dual problem 



(2.3) 



"itt - u lxx + Aui + Cu 2 = 
watt + u 2xxxx + Bui + Du 2 = 
Ul (t,0) = ui(t,%) = 
u 2 (t,0) = u 2 (t,n) = 
u 2xx (t,0) = u 2xx (t,ir) = 
ui(0, x) — uiq(x) and "it(0, x) = 



in M x (0,7r), 
in R x (0,7r), 
for t G R, 
for t G R, 
for < G E, 
Mil (a;) for x € (0, tt), 



u 2 (0,x) = M 2 o(a ; ) and u 2t (0,x) = u 2 i(x) for x G (0,7r), 

k V>(*) ="i (*,£)+ "2 (*,£) forteR. 

We prove that under some conditions this dual problem is observable. Feedbacks 
of this type are important for the engineering applications: as we will show in this 
paper, on various numerical aspect of these feddbacks, and we can see the works of 
Bourquin et al. [6] on physiscal experiences. 

3. Proof of the main results 



We consider the abstract observability problem (2.3), if the initial data are given 
by the formula 



and 



"10 



"20 



(x) = afc sin kx, 



fc=i 



(x) = ak sin kx, 



iti 



"21 



i (x) = bk sin kx, 



fc=i 



(x) = Pk sin kx 



fe=i 



k=l 



4 



ALIA BARHOUMI AND ABDELKADER SAIDI 



with only finitely many non vanishing coefficients a k ctk /3fc, then a simple 
computation shows that 



ii(t,x) = y2(c k e m +c_ fc e ikt )sinkx 



k=i 



and 



with 



and 



"2 



(t, x) = J](4e ifc2 * + d_ fc e _ifc2 ') sin kx 



k=l 



1/ .&fcx 1 b fc 

Cfc = - *y)> c_ fc = -(a fc +j — ), 



d * = 2^ fe fc2" d^ k = -{a k + i — ). 

If T > 2ir, then using Paxseval's equality and a result of Haraux [10] it follows 
that 

(Mt,oi 2 + k(f,77)| 2 )<ftx 



oo 

^ (K| 2 + fc- 2 |6 fe | 2 ) sin 2 fc£ + (|a fc | 2 + fc^ 4 |/3 fc | 2 ) sin 2 k V . 
fe=i 

It can be rewritten in the form 

(3.1) / T |«i(t,Ol a + MM)l a dt 

Jo 



We rewrite (2.3 1 as a first-order system 
(3.2) [/' = A*U, 17(0) = Z7 , V = 

by setting 

E/ := (wi,U 2 ,Ult,M 2 t); 

:= (itiO)W20)«ii!«2i), 
.A*(tti, U2, vi, Ua) : = (vi, «2, Aui — Ami — Cu 2 , — A 2 u 2 — Blti — DU2) 

and 

B*(u 1 ,u 2 ,v 1 ,v 2 ) := (ui(0, "2 (??))■ 

We introduce the dual space of the Hilbert spaces H^ )V denoted by H'^ := D° x 
£>° x D- 1 x L>- 2 , g := K 2 and we define the domain of definition of the linear 
operators A* and B* by 

D(A*) = D(B*) = D\ x L> 2 x .0° x £>°. 



Proposition 3.1. the system (2.3) verify the following four assumptions: 

(HI): The operator A* generates a strongly continuous group of automor- 
phisms e tA in H'^ . 

(H2): D(A*) C D(B*), and there exists a constant c such that \\B*Uo\\g < 
c\\A*Uo\\H>i, v forallU €D(A*). 
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(H3): There exist a non degenerate bounded interval I and a constant ci such 
that the solutions of (3.2) satisfy the inequality 



\\B*U\\ L 2( I; g) < C/||{7 ||h'£,77 

for all U a eD{A*). 
(H4): There exists a bounded interval I' and a positive number d such that 



the solutions of (3.2) satisfy the inequality 



\\Uo\\ n <c'\\B*U\\ LHl , g) 



for all U G D(A). 



Proposition 3.1 is an application of the abstract Komornik's method [15] . See 



also [5] for a proof of this proposition. 

Proposition 3.2. Assume (HI) — (H3) and fix T > arbitrarily. For any given 
Xq = (2/10,2/11,2/20,2/21) € and u = (vi,v 2 ) G L 2 (0, T; M 2 ), the problem ( fi~I] ) 

has a unique weak solution X — (2/1, 2/it, 2/2, 2/2t) G C(0,T, Wf )77 ), and the linear 
mapping (Xq, u) h- >• X is continuous with respect to these topologies. 

Proof. We rewrite ( |1.1[ ) as a linear evolution problem 



(3.3) X' = AX + Bu, X{0) = X , 

We can see [TS] for the necessity of the assyptions (HI) — (H3) in the abstract 
form. Next we define the solution of (1.1) by transposition. Fix Xq G 'H{, ?7 and 
u = (vi,V2) G L 2 (0,T;K 2 ) arbitrary. Multiply the equation (3.3) by the solution 
U of the equation in (3.2). Integrating by part formally between and T G K, we 



easily obtain the identity 
(3.4) 

<X(T),U(T) >u^,u> 



-< Xq. Uq >« t ,„,W^, ) 



+ / < u(s),B*U >j 



ds. 



Hence we define a solution of (3.3) as a continuous function X 



satisfying the identity (3.4) for all G Ti^ and for all T G R. This definition is 
justified by the following lemma. 

□ 

Lemma 3.3. As sume (HI) — (HA). For any given X , e H^ v andu £ L 2 (0,T;R 2 ), 
the problem (3.3) has a unique solution. Moreover, we have the estimates. 

( 3 - 5 ) ||^|U°°(0,T;M s ,„) < ^t(|!^o||« 5 ,^ + \\u\\ L 2(0,T;R2) 

with some constant Mt which does not depend on the particular choice of Xq and 
for all T > 0. 

In order to prove the stabilization estimate, we need to recall a general result 
proved in [T5]. Fix two numbers T > U) > 0, set T w = T + (2a;)" 1 , define 



and set 



eui(s) 



2uje 



-2u)T 



(T u 



if < s < T, 
s) ifT<s<T u , 



H'fli 



eul (s)(Be SJi Uo,Be sA U )g ds 



Then A w is a self-adjoint, positive definite isomorphism A u G L(H,H)'). Let us 
denote by J : Q — > Q' the canonical Riesz anti-isomorphism. 

The following result is a special case of a theorem obtained in [T5] . 
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Theorem 3.4. Assume (Hl)-(H^) and fix uj > arbitrarily. Then the problem 

(3.6) v' = {-A-BJB*AZ 1 )v, w(0)=«o, 

is well-posed in W. Furthermore, there exists a constant M such that the solutions 
of ( |3.6| satisfy the estimates 

(3.7) \\v{t)\\ n < M\\v \\ n e- Ut 
for all vq € "H and for all t > 0. 

In other words, this theorem asserts that the feedback law 



(3.8) 



W 



-JB*A: 



uniformly stabilizes the control problem 

v' = -Av + BW, v(0) = v 

with a decay rate at least equal to ui. 

The well-posedness means here that (3.6) has a unique solution v e C(M.;'H) for 
every v E W. 

3.1. proof of theorem |2.3[ Since hypothesis (HI) — (HA) are all satisfied, we 
may apply theorerr j3.4| In order to write down explicitly the stabilization result, 
we multiply the equatio n by u and we integrate by parts as follows (we use all 
conditions in JO} and (2.3)). This shows that if we write (2.3) in the form (3.3), 
then its dual (3.2) corresponds to (1.1). Furthermore, writing the operator 



A- 1 : (D*)' x (D-'Y x (£>«)' x (D~ )' 
in the matrix form 



D- 1 xD°x D- 2 x D° 



a; 



(- Pl 
A 2 i 

Aai 
\A 41 



Qi 
A 22 
A 32 
A 42 



Ai 3 A i4 \ 

-Pi Q 2 

A33 A 34 

a 43 A 44 y 



we have 

(3.9) ( Vl (t),v 2 (t)) 



-((Piyu + Qm)(t, 0, (P2V2t + Q2V2)(t, r,)). 



4. Numerical Approximation 



To perform a numerical computation we use a Faedo-Galerkin method. This 
allows us to approach numerically the operator A w using the family of the functions 
Wk(x) := \/2pK sin kx. An approximate solution y N — (y\ ,y 2 ) of the coupled 



system (1.1 ) is a solution of the problem 



(4.1) 



Vim 



V2,tt + V2,xxxx 



By? + Dy$ 
yf(t,0)=yf(i,7r)=0 
y 2 v (t,0)=yf(t,7r)=0 

y2,xx(t,o)^y2 N ,xx(t,^) = o 

V?(0,x) = y N {x) and |/{f t (0,i) 
1 1^(0, 1) = l$(a;) and < t (0,x) 



in (0,T) x (0,tt), 
in (0,T) x (0,tt), 
forte (0,T), 
forte (0,T), 
forte (0,T), 
for x e (0, 7i"), 
for x e (0, 7i") 
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In order to compute the state feedback law v N (t) — (v^(t),V2 (t)) — F (y N {t) , d t y N (t)) , 
an adjoint state is introduced : let u N (s 7 x) — (uf(s, x), {s 7 x)) be the solution 
of the coupled adjoint system : 



(4.2) 



L lss 
L 2ss 



,N 



Cu* = 



' u 2xxxx 



Buf 



<( S ,0) =<( S ,7t) 

<( S ,0)=<( S) tt) 



a 2xx 



(3,0) 



,N 



( S ,7T)=0 



Un-, (x) 



in (0,5) x (0,tt), 
in (0,5) x (0,tt), 
for s G (0,5), 
for s G (0,5), 
for s G (0,5), 
for x € (0, 7r), 
for x € (0, 7r),. 



(0,x) = Ui (x) and wf t (0,x) 
,v>2 (0,x) = u^o(x) and u^ t (0,x) - a 21 y 

Where s denotes a fictitious time , an d S a fictitious time horizon. The solution 
U N (s, x) = (u^(s, x),v,2 (s, x)) of (4.2 1 depends linearly on initial conditions U$ := 
(u^q, U20, wfi, u 2i)- Hence one can define the app roximate bilinear controllabiliy 
gramian for any solution U N — (u^ ) of ( |4.2[ ) with initial conditions 

(u 



r,N ~N ~N ~N 
10' u 20> ' 



xx,U2i) and parameter uj as : 



(4.3) 



N A N ) ■■= 



e- 2 " s «(s, (*, + v?(s, r,)v?(a, V )) ds. 



As pointed out in [6] we consider, from numerical point of vue, the function e~ 2ujs 
initially introduced in [2] and wich gives similar results as the general theory [T5] . 
For any Z = {zio, Z20, z\\, Z21} let Uq — {u^ , w^,, u^, u^} be the unique solution 
of the variational equation : 

(4.4) a» s (U?, Z) = ^ [(z u u& - z 10 u&) + (* 2 i<4 - zjjo«^ ))] dx. 

Jo 

Let know define the operator Cn by {u(1q, u^q, u^i, m^} = £^({^107^20^117^21}) 
and let also define : Pi : M. 4 — >• R and P2 : M 4 — > R respectively as the projection 
on the first and second component, i.e. Pi ({a, b, c, d}) — a and P2({a, b, c, d}) = b. 
The feedback law v N (t) — (v^ (t), v£? (t)) is then given by : 

(4-5) <(i) = -PiC N {y», y», (yf )', {y^)'}{x = 

and 

(4-6) <(f) = -P 2 £ N {y?, y?, (yf )', (y 2 v )'}(x = 77)) 

To compute the feedback vn (t) we use the expansion of initial conditions : 



iV 



A' 



l io( x ) =J2a° k sinkx, 



l n ( x ) = a k sm k x 



and 



k=l 



N 



k=l 



(x) = ^/^sinfcir, 



l 2i( x ) = ^2 Pi sinkx. 



k=l 



k=l 
Ni 



an approximate solution of problem (2.3) is then given by u^(t, x) — y\ —1 aj~(t) sin kx 



and U2 {t,x) = Ylk=i @ k (*) sm ^ x • With : 



and 



AW = ^(/3 fe -4 )e ' 



k2t + l(P°k+i§)e- 



■i k I 



ik z t 
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by setting rrii = sin(i£) and m = sin(irf) the second hand side of (4.3 1 is given by : 

N N S 

( 4J ) / e ~ 2U1S a i( s "> mim i + ^( s ) &OH™i) ds 

Let us know define {<P%(t),i^%(t),(l)}^{t),ilJf^(t)} as a solution at each time i > 0, 
of : 

(4.8) {0Sr(*),^W,^(i),^(*)} = ^{»fW,»f(*).(l^ (*))', (*))'} 



We can then write the operator £jv in a matrix form by setting : 

K^w&Nit), ^(*), ^ (*), = {(yf (*))', (% w (*))', -yf (t), -y 2 w (t)} 



Where the matrix K u ^n is defined by relation (4. 3), (4.7). Let us know define the 
inverse K~ l AN by a block matrix : 

K-^ N = (K^),a,p = l.A. 
the control law is then given by : 

N N 

v N (t) = (^(t),<w) = E(^w,5>oV*) 

ic=i fc=i 

Where 0^(4) = X n a'(t) - A' 12 a(i) and ^(t) = K 23 /3'(i) - X 24 /3(t) 

We compute the solution of the coupled system for N = 7 and value of parameters 

\/2 

A = 3. B = 2, C — 1 and D = ~. The damping points are chosen as £ = — and 

2 3 

i] = We show the efficiency of the feedback for two values of parameter lu in 

figure 1 (cj = 1) and figure 2 (uj = 10). 




figures 1 and 2 : solution of the coupled system , to = 1 (left) and w = 10 (right). 
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figure 3 and 4 : Energy of the coupled system , u> — 1 (left) and u> = 10 (right). 



In figures 3 and 4 we compare the energy of the coupled system in the space 
D® x x D® x D 2 ( solid line ) to the classical energy in the natural space 
£ 2 (0,tt) x H%(0,ir) x £ 2 (0,7r) x (i? 2 (0,7r) n H$(0,n)) (dashed dot line) for two 
values of the parameter uj = 1 and uj = 10. The computations shows that we have 
also uniform decay of the energy in the natural space. 
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